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>>: 

j^ ^ Abstract. A non self-similar change of coordinates provides improved 

matching asymptotics of the solutions of the fast diffusion equation for 
large times, compared to already known results, in the range for which 
1/^ , Barenblatt solutions have a finite second moment. The method is based 

on relative entropy estimates and a time-dependent change of variables 
_ which is determined by second moments, and not by the scaling corre- 

'"j , sponding to the self-similar Barenblatt solutions, as it is usually done. 



in 



X 



1. Introduction and main results 



Consider on M the fast diffusion equation 

O^ ' for some m G {rric, 1) with rric '■= {d — 2)/(i. Assume that the initial data 

is a given nonnegative function vq in L^iW^). It is weh known (see for 
instance [2j ) that the large time behavior of the solution is captured by the 

^ ■ Barenblatt solutions given for any (r, y) G R+ x M'^ by 

o' ^_ ^ 

B(T,y) := (1-Fr) m-,n, ^ Zj + 2d(m-mc) (^ + '^) '*^'"^'""' 1^1 



where Z) > is determined by the condition H't'ollLimd) = I|B(t, •)||2^i( 
Using entropy methods, it has been established in [5] how a linearized prob- 
5^ \ lem involving the relative entropy and the relative Fisher information de- 

termines the best rate of convergence towards the Barenblatt solution. The 
note [6] is devoted to a refinement of the estimates in which the depen- 
dence on t^o is clarified and the precise value of the best possible rate of 
convergence is computed in terms of a spectral gap of the linearized opera- 
tor associated to the relative entropy and the relative Fisher information for 
all values of m < 1. By taking advantage of the translation invariance, it is 
moreover possible to impose that the solution evolves in the orthogonal of 
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the eigenspace associated to the first non-zero eigenvalue of the linearized 
operator for m £ {mi, 1) with mi := (d — l)/d, thus providing an improved 
rate of convergence. The corresponding conserved quantity is the center 
of mass, while the generators of the eigenspace are the derivatives of the 
Barenblatt solution with respect to each of the coordinates. There is no 
other conserved quantity known, so that further improvements cannot be 
achieved directly by this method. 

One may however notice that the eigenspace corresponding to the second 
non-zero eigenvalue in the range m £ {m,i, 1) is generated by the infinites- 
imal dilation of the Barenblatt solution. It is therefore natural to try to 
adjust the Barenblatt solution by a scaling. This can be done by taking 
a time-dependent change of variables where the scale is determined by the 
solution itself and not anymore by its asymptotic, self-similar behavior, thus 
providing improved convergence rates. Asymptotically, we will recover the 
self-similar profile, but with a better matching. There is a price to pay: 
the rescaled equation has a time-dependent coefficient, which converges to 
a constant. From the point of view of the entropy - entropy production 
inequality, however, nothing is changed, which is the main observation of 
this paper. 

Let fhi := d/{d+2). In the range m, G (mi, mi), the infinitesimal dilation 
of the Barenblatt solution generates the eigenspace corresponding to the 
first non-zero eigenvalue. Our time-dependent change of variables therefore 
improves on the rate of convergence for any m G (mi, mi), and also for any 
m G (m-i, 1) if the center of mass is chosen at the origin. 

The reader interested in understanding the heuristics of our approach is 
invited to go directly to Section [2j In the remainder of this section, we will 
give a precise statement of our main result and some additional references. 

Define the mass and the center of mass of vq respectively by 

M := vody and xq := —- yvody. 

Consider the family of the Barenblatt profiles 

(2) S^(x) :=a-5 (C7M + i|xp)^ V x G M'^ 

where a is a positive real parameter and 

2(l-m) /d(m-mc)\ 

\ M* / _ md 



Notice that ||-Bo-||ii(iRd) = M for any cr > and B^j is a solution of 
V • \b ffji^'"^-™) VS"^-! - 2x11 = . 
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Let us recall the definition of rric-, mi, rhi and introduce the exponents m2 
and m2, for later use: 

d-2 d-l ^ d 

rric = — ; — , mi = — ; — , mi 



d ' ' d ' ' d+2 ' 

d+1 , ^ d + 4 

'^''■=dT^ ""^ '^''■=dT6^ 

which are such that m,c < mi < m,2 < 1, m-c < fhi < m2 and, if d > 2, 
fhi < mi and m2 < m,2- For later purpose, it is also convenient to define 

(3) Km:= [ \x\^Bidx= ^^~'^l^' MCM 

jRd m-mi 

for any m G (mi, 1). Notice that Km is finite if m > mi. 

If u is a solution of ([T]), consider the time-dependent scale R{t) defined 
by 

(ilogi? / Km \lim-m.) 



The justification of such a choice for r i— t- R{t) will be made clear in Sec- 
tionO Note that, in view of the fact that the quantity J^a \y — xo\'^ v{t, y) dy 
is non-decreasing in time along the solution of the fast diffusion equation ([1]) 
and because of its asymptotic behaviour, the time-dependent scale R{t) is 
such that log R{t) is increasing from zero to infinity. We also define o" as a 
function of r by the condition 



'^('^) •= ^^ — ^r\2 \y-xo\ v{T,y) dy . 
Km R{tY J^d 

As a consequence, the equation for R can be rewritten as 
(5) 2(j-i("'-™=) =i^^"'^(^~"') — 

and we can define for any x = y/R{t) G M the Barenblatt type solution B 

by 

^^'^''^^'^ WrY^V'Wr)) ^^^^^ 'S(r,x) :=B^(^)(2;) . 

The difference of B with the Barenblatt solution B is that B depends on 
j^d\y\^ vij^y) dy, so that they are only asymptotically equivalent, as we 
shall see later. The point is that B provides a better asymptotic matching 
than B. Our goal is indeed to measure the rate of convergence of u towards B. 
For this purpose, it is convenient to change variables and study the rate of 
convergence of u defined by 

u(r, x) = R{tY V (r, xq + R{t) x) 
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towards 55. Let us consider the relative entropy of J. Ralston and W.I. 
Newmann defined in [221 [24] by 

£{t) := -^— f [n™ - QS"^ - m ^""'^ {u - 55)1 dx . 

Theorem 1. Assume that m G (mi,l), d > 2. Lei v he a solution of ([T]) 
wit/i initial datum vq G LJ',_(]R'^) such that v^ and y i— )• \y\'^ VQ{y) are inte- 
grable. With the above notations, we have 

limsupi?(r)^('").S(r) < oo , 

T— >00 

where R{t) 



J- d [m — mc) 


as r — 7- oo and 






r ((d-2)m-(d-4))2 
4(l-m) 


i/m G (mi,?7i2 


y{m) = < 


4:{d + 2)m-id 


i/m G [m2, w,2 




. 4 


i/m G [m2, 1) . 



Once a relative entropy estimate is known, it is possible to control the 
decay rate of n — *B in various norms, for instance in L'^(M , dx) for 



^ fi 2d{l-m) 1 



or in C'^, by interpolation. Up to a change of variables, this also allows to 
prove decay rates of v — B. See [5] for more details. 

Compared to the results of [6], an improvement for any m > fhi has been 
obtained. The values obtained in [6] for 7(771) are indeed 

{2 dm — 2 {d — 2) if 771 G (mi,77ii] , 
2 if 777 G [7771, 1) , 

except that in [6] the scale R{t) is determined by the self-similar Barenblatt 
solutions (both scales are anyway equivalent as r — )• 00: see Lemma[8]). Also 
see Figure 2 at the end of this paper for more details on 777 i—)- 7(777) in the 
setting of [6j compared to the results of Theorem [TJ 

Compared to other methods, it may look surprising that the scale R{t) 
and, as a consequence, the coefficient a both depend on the solution t; of ([T|). 
Asymptotically, as r — )• 00, R{t) is equivalent to the scale given by the self- 
similar change of variables, but what has been gained is a better matching 
with the closest Barenblatt solution. The family of the Barenblatt solutions 
is globally invariant under scaling and, among all such solutions, there is 
one which is closer to our solution of the evolution equation: the one with 
the same second moment. 

Convergence results of a suitably rescaled flow associated to ([T]) towards 
an asymptotic profile has been established in [19] for 777 > 777c (also see 
for instance [26]) and in [U [15] for m < rUc- Getting rates of convergence 
beyond a simple interpolation between mass and uniform estimates has re- 
quired the use of the relative entropies introduced by J. Ralston and W.I. 
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Newman in [22^123] . First results in this direction have been achieved in |13j 
using the entropy / entropy-production method of D. Bakry and M. Emery 
(also see |llj for general diffusions and [Ij for an overview) and in [16] us- 
ing sharp Gagliardo-Nirenberg interpolation inequalities. F. Otto made the 
link with gradient flows with respect to the Wasserstein distance in [23], 
and D. Cordero-Erausquin, B. Nazaret and C. Villani gave a proof of the 
corresponding Gagliardo-Nirenberg inequalities using mass transportation 
techniques in p^. The condition m > mi was a strong limitation to these 
first approaches. Gagliardo-Nirenberg inequalities indeed degenerate into a 
critical Sobolev inequality for m = mi, while the displacement convexity 
condition also requires ttt. > ttt-i. For m < mi, various limitations appear. 
To work with Wasserstein's distance, it is crucial to have second moments 
bounded, which amounts to request m > rhi for the Barenblatt profiles; see 
for instance [17^ [T8] . Linearization of entropy estimates around the Baren- 
blatt profiles has been considered in [12^ [20j for ttt, G (mc,l). In a certain 
sense, this is also the strategy in [T71 [18] . Integrability of the Barenblatt 
profiles means m, > m,c- This condition has been removed in a series of 
recent papers (see [H [3 El [7] ) together with a clarification of the strategy 
of linearization of the relative entropies, at least from the point of view of 
functional inequalities. In this paper, we shall however restrict ttt. to the in- 
terval {fhi, 1), for spectral reasons that are explained in Section[3|and for the 
second moment to be well defined. For m < mi, even with an appropriate 
definition a relative second moment, our method gives no improvement on 
the convergence rates because of the presence of the continuous spectrum. 

Rescalings and convergence towards Barenblatt solutions, or intermediate 
asym,ptotics, has not been the only issue of large time asymptotics. We can 
for instance quote [25] for a study (in the porous media case) of the time 
evolution of the second moment, and [1011181 [211 [9] for the search of improved 
convergence rates when moment conditions are imposed in the framework 
of Wasserstein's or other Fourier based distances. The question of improved 
rates has been precisely formulated in [18], and solved in [6] in the weighted 
L^ framework that we shall use in this paper, as far as the first moment 
(position of the center of mass) is concerned. The main contribution of this 
work is to explain how improvements based on the second moment can also 
be achieved. 

This paper is organized as follows. In Section [21 we explain how faster 
convergence results can be achieved by introducing an appropriate time- 
dependent rescaling, which is given by ^ and not by the explicit depen- 
dence of the Barenblatt solutions. Improved Hardy-Poincare inequalities 
are established in Section [3l using the spectral results of [17, [ISj and the 
spectral equivalence found in [6] . The large time behaviour of the solution is 
studied in Section HI The proof of Theorem [1] is then completed in Section [5j 
Further considerations on the case d = 1 and the limiting regime as ?n, — )• 1_ 
are presented in the last section. 



6 j. dolbeault and g. toscani 

2. The relative entropy approach 

The result of Theorem [1] is easy to understand using a time-dependent 
rescahng and the relative entropy formalism. Define the function u such 
that 

(6) v{T,y + xo) = R-'^u{t,x) , R = R{t) , t = ilogi?, x = ^ 

R 

where u is a solution of ([T]) with initial datum vq G L^(M ). A simple 
computation shows that u has to be a solution of 



,„. du „ 



^f(m-m,)y^m-l_2 3. 



t >o, xe 



with initial datum uq = vq (we assume that R is chosen such that -R(O) = 1) 
and a given by 

2|^-f(m-nic) _ pl-d(l-m) "-R 

dr 

which is nothing else than ^. By virtue of the definition of R{t), the new 
time t = 2 log R{t) increases monotonically from to +00. Consequently, 
the old and new times can be uniquely related, and r can be expressed in 
terms of t through the inverse function of R, so that r = i?~^(e^*). Using 
this transformation and with a slight abuse of notations, we shall consider 
from now on o" as a function of t. It is important to notice that, as long as 
-^ 7^ 0, the Barenblatt profile B^j is not a solution of ([7]), but we may still 
consider the relative entropy 

J^M ■■= -^ I [u'^-B^-m B^-' {u - B,)] dx . 

Let us briefly sketch the strategy of our method before giving all details. 

If we consider a solution of ([7]) and compute the time derivative of the 
relative entropy, we find that 
(8) 

l--)i"<'- ■)! - ^ [i ->i) ,„,.r ^ L (-•"- - -:i.r ) t - • 

Here comes the main difference with previous works. As we shall see below 
in the proof of Lemma [2] (also see Remark [1]), the first term of the right 
hand side in ([8|) involves 

i^J^.H] =!!12^([ \,\^B^dx- [ \x\'udx) 

where D{a) := A((j'^(»"c-m)/2^^ When taking a time-dependent rescaling 
based on the self-similar variables, one finds that a is constant in t, so that 
^ = and the term -^ -^ {Ta[u\) does not contribute. In our approach, 
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a depends on t but can still be chosen so that this term does not show up 
either. It is indeed enough to require that 

(9) / \x\'^ Bcr dx = / \x\'^ u dx , 



which amounts to ask that R solves the ordinary differential equation (j3]), 
to obtain that ^ (Jv[u]) = 0. This will be justified in the first step of our 
method, below (see Lemma [2]). 

In a second step, we shall use the fact that 

d. 



(10) -^.(,)[n(t,.)] = T^T^^/ - 



V 






2 



dx . 



From there on, the computation goes essentially as in |5J[6]. For complete- 
ness, we will briefly reproduce it. However, with our choice of a, we gain 
an additional orthogonality condition which will be explicitly stated in the 
third step of the method: see Lemma El This orthogonality condition is 
the crucial point (see Corollary [7]) for improving the rates in Theorem (H 
compared to the results of [6J. Now let us give further details. 

First step: choice of the scaHng parameter. For a while, we do not 
need to take into account the dependence of a in t. The main idea of this 
paper is indeed to choose a in terms of u by minimizing a i— t- Jv [u] , so that 

d ^ , , 
da 

Lemma 2. For any given u G -L^(R'^) such that n™ and |xpu are both 
integrahle, if m £ (mi,l), there is a unique a = a* > which minimizes 
a I—)- J>[n], and it is explicitly given by 

(11) a* = -^ f \x\'^udx . 

For a = a* , the Barenblatt profile B„* satisfies ([9]). 

In (jlip . Km is the constant which has been defined in ^. The condition 
m > fhi guarantees that B™ is integrable and Km = f^d \x\'^ Bi dx is finite. 

Proof. We have to minimize 

e(a) := (1 - m) / B^ dx + m [ B^'^ u dx , 

which corresponds to the two cj-dependent terms in the expression of Jv[u]. 
Using the fact that ^^-^(x) = ct'^(i-™)/2(Ca/ + \x\'^/(t) and J^, udx = M, 
it is easy to see that 

e{a) = (1 - m) al ^^-""^ [M Cm + Km] 



|2 

+ ma^'^'''^> [MCm+ I '-^udx 



±{l^m) ( I \^\ 



a 
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Collecting terms, we get 

e(a) =a5(i-'") [MC7M + (l-m)i^M] +"1^1(1-"^)-! / \x\'^ u dx 

Hence, by optimizing on o" > 0, we find that a = a* is given by 

m{m-mc) /" i |2 , 

(^ = -, TT T^ -, ^ t" / UE U ax 

{l-m)[MCM + {l-m)KM] 
Using ([3]), we observe that 

m {m — rric) 



(1 - m) [M Cm + (1 - m) Km] Km ' 

which proves Q and therefore pT]) . D 

Remark 1. To prove ([9]) directly, we may notice that a* is determined by 
the condition 

= e'{a) = m (1 - m) f B^-^ ^ {B^ - u) dx . 
J^d da 

Hence we obtain 

= m / (B^-u) -^ — dx = m {B^ - u) [C{a) Cm + D{a) jxH dx , 
jRd do- J^d \ J 

where C{a) = i|(fT'^^^"™^/^), and D{a) = £{o'^^^~"'^^^'^) ■ Taking into 
account that both Bu and u have the same mass M , we get ([9]). 

The dependence of a in i when u is a solution of d?]) has not been taken 
into account yet. The choice of Lemma [2] determines an ordinary differential 
equation for R in terms of Jg^j \x\^u{t,x) dx. Undoing the time-dependent 
rescaling ([6]), this equation is exactly dH. With the choice R{Q) = 1, we 
recall that 110 = ^0- 

As already mentioned, the choice of a in Lemma[2]has a major interest. If 
we consider a solution of ([7]) and compute the time derivative of the relative 
entropy, we find that the first term of the right hand side in ([8]) drops so 
that 

and we are back to the usual computations in self-similar variables. 

Second step: the entropy / entropy production estimate. Accord- 
ing to the definition of Sg-, we know that 

2x = cji('"-'^^)VS™-^ 
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Using ([7]), we obtain 

,m—l 



m 



m — 1 



B-^-A ^ dx 



ma{t) 



^(m-mc) 



V 



,m— 1 



B 



»n— 1 
ait) 



dx , 



1 — m 

thus proving (jlOp . Let w := u/B„ and observe that the relative entropy can 
be written as 



-7>N 



m 



w 



1 — m 
Define the relative Fisher information by 



1 (w" 

m 



1) 



1a[u\ :-- 



1 



m 



-V[{w' 



m— 1 



i)b: 



771—11 



B':!}dx . 



BfjW dx 



For a sohition u of ([7]), we find that 

I J-^(,) [t.(i, .)] = - m (1 - m) a(t)i(— ^) X,(,) [^(t, •)] V i > . 

As in [5| Lemma 3] (also see ^) we can estimate from below and above the 
entropy Fa[u] by 



2 D2-m 



dx 



(12) /i™-2 / |/|2S2-- dx < - J-^M < h-"'^ f \f\'Bl 

where f := {w — l)B^~^, hi(t) := inf]gdit;(t, •), h2{t) := sup^dw{t,-) and 
h := max{/i2, 1/hi}. The fact that /i is bounded for any f > is easy to 
prove by the Maximum Principle if h{0) is finite. See for instance [5] for 
more details. Even if /i(0) is infinite, h is anyway bounded for any i > 0, 
large enough, when m G (rn,c, 1); see [HI Theorem 1.2]. By [5, Corollary 1], 
we also know that lim^-^oo h{t) = 1. 

According to [5< Lemma 7] (also see [6]) , the generalized Fisher informa- 
tion satisfies the bounds 

(13) / \Vf\^B^dx<[l + X{h)]lM+Y{h) [ \f\^Bl-^dx 
where 

d{l-m) \C0 (2-771) _ J <d{l-m) [h^ (2-m) _ 1] =. Y{h) . 
Notice that X{1) = Y{1) = 0. 

Third step: orthogonality conditions. To obtain decay rates of i i— t- 
^a{i) [^(^; ■)] is now reduced to establish a relation between J^^ |/p ij2-77i ^^ 
and Jj^d I V/p i?o- dx. This is the purpose of the next section, but before let 
us make a few additional observations on the properties of / = i?™~2 ^ _ 
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Lemma 3. Let u be a solution of ([7]) and f = B^^^ {u/B^ — 1) where 
a = a{t) is defined by ([5]). With these notations, the function f has the 
following properties, for any t > : 

(i) Mass conservation: f^a f{t,x)B'^~"^ dx = if m > rUc- 
(ii) Position of the center of mass: J^dX f{t,x) B'^~"^ dx = if m > 

(d-l)/id + l). 
(iii) Conservation of the second moment: J^^ |xp f{t,x) B^ "^ dx = if 
m > nil. 

Notice that Property (ii) has already been used in [6l Theorem 7] to 
obtain improved rates of convergence for m £ (mi,l). Property (iii) is 
new and arises from the fact that the change of variables ([6]) is chosen in 
Lemma [2] by imposing a condition on the moment and not according to the 
self-similar variables corresponding to the Barenblatt solutions. Restrictions 
on m are such that mass, first and second moments are well defined for B^j. 
As in [^ , such conditions could be relaxed by considering moment conditions 
on / only, or relative moment conditions; we will however not pursue in this 
direction, as such an approach does not improve on the rates of convergence. 

Proof. It is straightforward to rewrite the conservation of mass along evolu- 
tion into Property (i). 

Equation ([1]) being independent of y, it is clear that Jj^^ y w(t, y) dy = 
J^dyvo{y) dy for any r > 0. As a consequence, the center of mass of u is 
located at x = 0: J^d xu{t,x) dx = 0, and so we get Property (ii). 

Finally, Property (iii) is a direct consequence of Q. D 

3. Improved Hardy-Poincare inequalities 

When M = M* and a = 1, with the notations of Section [H the quantities 
Jjgd |/P-B^~"^ ^^ ^^^^ /iRd |V/pi?CT dx involve various powers of (1 + |xp). 
On M , we shall therefore consider the measure d^a '■= p-adx, where the 
weight /ia is defined by ^ia{x) := (1 + |xp)", with a = l/{m — 1) < 0, and 
study the operator 

J~-a,d ■= -fJ-l-a div [HaV-] 

on the weighted space L^ (dfia-i). The operator Ca,d is such that 



/ f{C^,df)dfi^_i= [ \Vf\^dfic 



Notice that in the range m G (?7ic,l), that is for a £ {—oo,—d/2), dfia 
is a bounded positive measure. If additionally m > mi, that is for a < 
— (d + 2)/2, then dfia-i is a bounded measure and J^d |a^P d^a is also finite. 

Let a* := -{d-2)/2. Based on [13 ig HE] , we have the following result. 

Proposition 4. The bottom of the continuous spectrum of the operator 
Ca,d on L? {dfia~i) is A™JJ* := (q — q*)^. Moreover, Ca,d has some discrete 



ASYMPTOTICS OF FAST DIFFUSION EQUATIONS 11 

spectrum only for a < a^. For d > 2, the discrete spectrum is made of the 
eigenvalues 

X^f^ = -2a{£ + 2k)-4k{k + £ + ^-l) 

with I, k = {),!, ... provided {I, k) / (0,0) andl + 2k-l < -{d + 2a)/2. If 
d = 1, the discrete spectrum is made of the eigenvalues X^ = k (1 — 2 a — k) 
withk GNn [1,1/2 -a]. 

Let a* := —{d — 2)/2. The following result has been established in [5]. 

Corollary 5 (Sharp Hardy- Poincare inequalities). Let d > 2. For any 

a G (— oo,a*) U (a*,0), there is a positive constant Aa^d such that 

(14) K,d [ \f\^dfla^l< [ Nfl^dfla yfeL^dfla-l) 

under the additional condition f^^ f dfia-i = if a < a^^. Moreover, for 
d> 3, the sharp constant A^^d is given by 



A 



' l{d-2 + 2af i/aG [-^,a*)u(a*,0) 
-4a-2d ifae [-rf,-^) , 

— 2a if a ^ {—CO,— d) . 



a,d 



For d = 2, Inequality (J14l) holds for all a < 0, with the corresponding 
values of the best constant Aq,^2 = ck^ for a E [—2, 0) and Aa.2 = —2 a for 
a £ (— oo,— 2). For d = 1, (fT4]) holds, but the values of Aa^ are given by 
Aq,i = -2aifa< -1/2 and Aa,i = {a-l/2f if a £ [-1/2,0). 

The constant A^^ ^^ is determined by the spectral gap and corresponds 
either to the lowest positive eigenvalue, Ai^o or Ao,i, or to the bottom of the 
continuous spectrum, A™JJ* := \{d + 2a — 2)^ (see Fig. 1). 

With additional orthogonality conditions, one improves on the spectral 
gap in the range for which discrete spectrum exists. A first result in this 
direction has been achieved in [6j for solutions with center of mass at the 
origin. Here we give a refined version of it, by going to the next order, that 
is, by considering functions with zero moments up to order two. 

Corollary 6 (Improved Hardy-Poincare inequalities). Under the assump- 
tions of Corollary\^ if a < —{d + 2)/2 and f £ L'^ {d^a-i) is such that 

f dfia-i = , / X f dfia-i = and / |x| f dfj,a-i = , 



then (|14p holds for any d> 3 with 

f i(d-2 + 2a)2 ^faG[-^,-^) , 
-8a-4{d + 2) ifae[-{d + 2),-!^] , 
-4a ifa£{-oo,-{d + 2)]. 

For d = 2, Aa^2 = o? for a E [—4, —2) and A^^i = —4a for a E (— oo, —4). 



A. 
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The constant A^^^ is now determined either by the lowest of the two eigen- 
values, Ai^i or A2,o, or by the bottom of the continuous spectrum, A™JJ* (see 
Fig. 1), since the components corresponding to the eigenspaces associated 
to Ai and Aq i have been removed. 




Spectrum of 
l-m)£i/(,„-i),,( 



Essential spectrimi 
of(l - m)£i/(„,_i)j 




Figure 1. Spectrum of Ca,d as a function of a (left), and 
spectrum of (1 — m) 'Ci/(^_i) ,^ as a function of m (right), for 
d = b. 

A crucial observation is that we can scale Inequality ()14p . 

Corollary 7. Let M > and d > 2. With B„ defined by ^, let f be a 
function in L^{B'^~'^dx) such that 

[ (l,x,|xp)/5^-™dx = (0,0,0) and Vf eL^{B^dx) . 

Then the inequality 

holds with Aa.d ols in Corollary {^ 
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Proof. The proof relies on a simple change of variables. From Corollary [Gj 
we know that 



A,,,/ Ig^l + \x\Y-' dx < f \Vgmi + \x\ 



2\o 



dx 



for any g ^ L^ (dfia-i) satisfying the conditions of Corollary [H Then Corol- 
lary [7] holds for / such that f{x) = g{x/\/C]\j a ) for any x G M*^, which 
concludes the proof. D 

4. Estimates on the second moment 

Up to now, we have not determined the behavior of R{t) as r — )• cxd, nor 
the fact that a has a finite, positive limit as i — ?> oo. These properties can 
be deduced for instance from [5]. For the convenience of the reader, let us 
give some details. As in [5], consider the standard change of variables 

(r,y) ^ (t = I logRo{T),x = ^^ ^ 

based on the self-similar behavior of the Barenblatt solution B(r, y), where 

1 
Rq{t) := {l + 2d{m — mc) r) <*('"-'"<:) . 

With the above change of variables, if f is a solution of ([T]), then the function 
U(t,x) ■.= v{T,y) solves 

BU 

— + V- [[/ (VC/"^-^-2x)] =0 t>0, xGM'^, 

which is nothing else than ([7]), except that here a is replaced by 1. It has 
been established in [5] that 

lim e^^^Fi[U{t,-)\ = 

t— >oo 



for some positive constant c. By Holder's inequality, we know that 

<( [ UB-^-^ dx] ( [ ff^dx 

KJR'i J KJR'i 

thus proving that 

J^i[U] > I B'^dx + ^^ I UBV"-^ dx 

jRi l-m J^d 

- -^ ( f UB"^-' dxY ( f BT dx) 
l-m \J^d ) \J^d J 

With 

^ f^.UB'^-Ux ^ f^,\x\''Udx + MCM 
k'iBTdx Km + MCm 



l—m 



1—m 
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this can be rewritten as 



hence showing 



l-m(k-l) 
m — 1 



~ |k(t) — 1| as t — )• oo , 



hm e'^*|k(t) - 1| = . 



As a consequence, we observe that 

/ \x\ U{t,x)dx — KM 



0{e ^*) ast—'oo, 



Undoing the change of variables, we find that 



1 



\y\'^v{T,y) dy- Km 



which, by definition of a, gives 
1 



0{Mt) 



-c/2^ 



as r ^- oo 



a 



(1 + 2 d {m - rUc) t) 



d {m — mc) 



1 + 0{Ro{t)-^/' 



R{tY 

as T — >■ oo, where r i— ;■ R{t) is given by (JH. Using ([5]), this means that 

1 1 



1 dR „ ,„ , ^2 _^-|(m-mc) 



R{t) dT 



2{Ro{rya) 



d [m — nic) T 



:i + 0(r-)) 



with e = niin{l,c/2} > 0. With these estimates, we can prove the following 
result. 

Lemma 8. With the notations of Sections\^and\^ 

1 
R{t) ~ r<*(™-'"c) ag 7" _). oo 

and, as a function oft, 1 1— )• a{t) is positive, decreasing, with 

lim a{t) =: a^o > . 

More precisely, we know that for some Coo > 0, 

1 

i?(T) = CooT''(™-™=) (1 + o(l)) as T -)- OO 
2 

and (Too = {2d{m — nic))''^"'^"^''^ /C'^. However, the value of cJoo in terms 
of vq is not known. 

Proof. The asymptotic behaviors of R and a are direct consequences of 
the above computations. We only have to prove the monotonicity of a. 
According to Lemma El we know that 



a{t) 



Km 



|a;p u{t,x) dx 
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Using ([7]) and integrating by parts, we compute 
\x\ u{t,x) dx 



da 1 d I 2 



dt Km dt 

ttiKm 
On the other hand, by Lemma O we have 

B^^^l u{t, x)dx= f S^,) dx 

SO that 

-^.(t)Ki,-)]= / f^dx 

./rod m — I 



^"^ i^^,^i(i— )-_^ / u^dx. 



d{l — m)^ 1 — m 

Hence we have proved that 

dt m Km 

which completes the proof. D 

Remark 2. We may notice that the second moment converges to a well de- 
fined value: 

This is not the one which is usually found by considering the time- dependent 
resettling corresponding to the self-similar Barenblatt solutions. 

5. Proof of Theorem [T] 

We are now ready to resume with the relative entropy estimates and 
conclude the proof of Theorem [TJ Using p^ , (fT3|) and Corollary [71 we find 
as in [6] that 

(15) TAu] < ^jp^±^mai^—^^XAu] 

2 [Aa,d-cry{h)\ 

as soon as < /i < /i* := min{/i > : Aa^d — supjg]jj+ (T(t) y(/i) > 0}. 
Two differences with [6j arise: A^^d has been improved in Corollary [71 to the 
price of a factor a, which however plays no role because it also appears in 
the computation of ■^F^(i-^[u{t, •)]. The m factor is present because Jv has 
not been normalized as in [6j , and also because the equation for v is not the 
same. 

As in [6], uniform relative estimates hold, according to [TJ Inequality 
(5.33)]: for some positive constant C, we have 



Q<h{t) -I <CF^i^t)Ht, ■)]"+'- 



{d+l)m 
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Suniniarizing, we end up with a system of nonlinear differential inequalities, 
with h as above and, at least for some t* > large enough. 






;i-m) j;(t)[n(t,-)] 



for any t > t^. Gronwall type estimates then show that 

limsup e2(^"™)^i/{— i),d*jr^(j^[u(t,.)] < oo 



4— >CXD 



Notice as in [5j that for some constant c > 0, limi_j.oo e'^ {h{t) — 1) = 0, so 
that the fact that the quotient Xo-[ii]/J-'o-[u] in ([T5|) is not estimated exactly 
by (2 AQ,^rf)/(m(T) plays no role for the rate of convergence. This completes 
the proof of Theorem [1] with 7(771) = (1 — tti) Ai/(m_i) ,^ (see Fig. 2). D 

Remark 3. Exactly as in ^ Corollary 3], explicit estimates of the constants 
can be obtained. If vq/Bi is uniformly bounded from above and from below 
by two positive constants, an estimate of \\TQ.sVip^^^R{T)"'^'^> Tfj(t)[u] in 
Theorem [7] can be given in terms of vq by computing a Gronwall estimate. 




Figure 2. For d = 5, the value of m 1— t- j{m) is given by the 
curve of Case 1 when no assumptions are made on the initial 
data. The curve of Case 2 corresponds to the case where the 
center of mass is chosen at the origin, as in [6], while the 
curve of Case 3 corresponds to the exponent found in Theo- 
rem [TJ 



asymptotics of fast diffusion equations 17 

6. Concluding remarks 

The case of the heat equation, i.e. m = 1, is not covered in our approach. 
However, we may pass to the hmit as m — )■ 1_ in CoroharylTl for the special 
choice l/o" = 2 (1— m). Both weights B'^~"^ and B^ converge to the Gaussian 
weight, so that the conditions of Corollary [7] become 

/ (l,x,|rE|2)/ e-l^l'/2(ix = (0,0,0) and V f € L^e'^^^' /^ dx) . 

By requiring the orthogonality with respect to all Hermite polynomials up 
to order two, we achieve the improved Poincare inequality 

2 / I/I' e-l^l'/2 dx< [ |V/|2 e-l-'l'/2 dx . 

Compared with the results of [6j, nothing is gained, as 

lim (1 — m) Aoi = lim (1 — m) A20 = 4 , 

m— >1_ 771— J-l- 

where \ek are defined in Proposition H] and a = l/{m — 1). See |3J for more 
details on improved convergence rates of relative entropies in case m = 1. 

For completeness, let us extend our results to the case d = 1, which is 
very simple. Eigenvalues of Ca.i are ordered uniformly with respect to q, 
according to Proposition HI Let a < —1/2 and consider f £ L'^ (dfia-i) such 
that 



f dfia-i = , / X f dfj,a-i = and / |x| f dfj,a-i = 

: Jr Jr 

Then ^ holds with 

r (a-i)2 ifae [-f,-^) , 



A, 



a,d 



6(a + l)) ifaG (-oo,-|]. 



Theorem 9. Assume that d = 1, m £ (1/3,1). Let v be a solution of ([T]) 
with initial datum vq G L^(M) such that u™ and \y\'^ vq are integrable. With 

the above notations, we have lim.snp^_^^ R{Ty^''^> £{t) < 00, where R{t) ~ 

1 
r"i+i as r —7- 00 and 



7(m) 



4(i-m) y "'' ^ U' 5J ' 
6m i/?n, E [|, 1) , 
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